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Рассматривается проблема потери устойчивости стационарных траекторий движения инструмента относи- 
тельно заготовки с учётом периодического изменения параметров взаимодействующих подсистем и пара- 
метров динамической СВЯЗИ, формируемой процессом резания. Изменения параметров обусловлено ПеЕрио- 
дическими вариациями жёсткости заготовки в направлении, нормальном к оси её вращения, а также вариа- 
циями припуска, периодическими с периодом её вращения. Показано, что периодические изменения пара- 
метров вызывают не анализируемые ранее механизмы потери уст ОЙЧИвоОсти, Причём, потеря ус ТОЙЧИВОСТИ В 
этом случае принципиально зависит от частоты вращения шпинделя. Кроме этого, потеря ус тоЙчивости свя- 
зана также с кинематическими возмущениями, обусловленными периодической составляющей скорости про- 
дольной подачи суппорта. Приводятся примеры потери устойчивости в зависимости от параметров матема- 
тических моделей системы резания. 
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Введение. Одним из направлений повышения точности обработки на металлорежущих станках 
является управление процессом на основе изменения траекторий исполнительных элементов 
станка [1]. При этом используется системный анализ и синтез [1, 2], основанный на использова- 
нии асимптотического поведения системы, динамика которой описывается нелинейными диффе- 
ренциальными уравнениями, имеющими малые параметры при старших производных. Таким об- 
разом, система состоит из двух иерархически связанных подсистем «медленных» и «быстрых» 
движений. Ранее рассмотрены проблемы асимптотической устойчивости траекторий в этих подси- 
стемах [3]. При этом не учитываются изменения свойств подсистем, связанные с периодическими 
изменениями их параметров в математических моделях динамической связи, формируемой про- 
цессом резания. Однако, периодические изменения параметров, период повторения которых 
определяется частотой вращения шпинделя, естественны в рассматриваемой системе. Они могут 
вызывать параметрические явления, приводящие к потере устойчивости траекторий движения 
инструмента относительно заготовки. Именно проблеме потери устойчивости системы за счёт па- 
раметрического возбуждения стационарных траекторий посвящена настоящая статья. 

Обоснование математической модели. Периодические изменения параметров линеаризован- 
ного в окрестности стационарной траектории уравнения динамики имеют две причины. Первая 
обусловлена периодическим изменением жёсткости в подсистеме заготовки, закреплённой в па- 
троне. В этом случае, как показывают экспериментальные исследования, значения жёсткости 
подсистемы заготовки в плоскости, нормальной к оси её вращения, периодически изменяются по 
углу вращения заготовки (рис. 1). В статье мы ограничимся случаем, когда деталь закреплена в 
трёхкулачковом патроне. Периодические изменения жёсткости обусловлены тем, что при распо- 
ложении заготовки против кулачка патрона жёсткость всегда больше, чем при её расположении 
между кулачками. При этом параметры жёсткости подсистемы заготовки, закреплённой в шпин- 
деле в направлении, проиллюстрированном на рис. 1, и в ортогональном к нему направлении, 
имеют примерно неизменные распределения, но они сдвинуты по фазе для трёхкулачкового па- 


2 ыы 
трона на угол 50 Характерно, что уровень модуляции параметров жёсткости подсистемы обра- 


батываемой заготовки зависит от геометрии кулачков и способа закрепления детали в заднем 
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вращающемся центре. Уровень модуляции определяется параметром и в выражении 
с(1+нс0$ 30). Кроме этого, частота @ — есть частота вращения шпинделя. Ранее показано 
[2], что параметры жёсткости подсистемы инструмента в двух ортогональных направлениях яв- 
ляются симметричными, поэтому матрицы жёсткости подсистемы заготовки можно представить в 
виде: 

с(1+нс0539) 0 


с= , 1 
0 стерео [ 30+ За] (1) 


так как при симметричной геометрии и условиях закрепления всякая ортогональная система ко- 
ординат является главной. Здесь жёсткость подсистемы заготовки рассматривается в плоскости, 
нормальной к оси вращения шпинделя, так как жёсткость в направлении оси вращения шпинделя 
обычно на порядок больше. 





Рис. 1. Пример изменения жёсткости заготовки, закреплённой в трёхкулачковом патроне, при изменении угла 
поворота шпинделя 


Вторая причина обусловлена погрешностями профиля поперечного сечения заготовки, 
вызванными неточностью установки заготовки в зажимном приспособлении и смещением оси за- 
готовки и оси её вращения. Кроме этого, на периодические изменения припуска влияют радиаль- 
ные биения шпинделя, которые, как показано в работе [4], также имеют круговые траектории в 
вариациях относительно идеальной оси вращения шпинделя. Все эти погрешности, представимые 
в виде ограниченного ряда Фурье, вызывают изменение элементов матрицы жёсткости процесса 
резания, регулярные по периоду вращения заготовки. Если полагать, что величина подачи на 
оборот в процессе резания остаётся неизменной, то, следуя работе [3], формирование двух со- 
ставляющих сил резания, зависящих от деформационных смещений инструмента в направлении 
подачи, определяется выражением: 





1=М 
Ед (ЕТ, я = ХС. | - Ум со | (27 - 1) ах; 
аР, = 
Рд (В) +Т, де = Ха | - У, с0$ | (2/- 009 4х; (2) 
аР, < 
д (ВТ, = Хо | +, соз | (27 — о ах, 
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где Ал (Е), Бл (Е), Руд (Е) — составляющие сил, обусловленные деформационными смещения- 
ми инструмента в направлении подачи, действующие соответственно в направлении, нор- 
мальном к образующей формируемой поверхности, в направлении скорости резания и скоро- 
сти подачи; х,, Х,›, Х. — угловые коэффициенты ориентации силы резания в пространстве; 
с, =рЁ, — приведённая жёсткость процесса резания, рассматриваемая без учёта вариаций 
припуска (здесь Ё,, р — глубина резания и давление стружки на переднюю поверхность ин- 
струмента); Т,, Т,, Г, — постоянные времени, определяющие запаздывание сил при измене- 
ниях ДАХ = Х, (Е) -Х. (Е-Т). Здесь Х, (Е) — упругое деформационное смещение вершины ин- 


струмента в направлении подачи. Деформационное смещение заготовки в направлении пода- 
чи не рассматривается по указанным выше причинам. 

Выражения (2) показывают зависимость формируемых сил резания от деформационных 
смещений с учётом вариаций припуска, периодических по периоду вращения заготовки. В частно- 
сти, если рассматривается абсолютно круглое тело заготовки, ось которого смещена относитель- 
но оси вращения шпинделя, то 


[= 
Е (ЕЁ) +Т, бл = Х,сь {1 + с0$ [9Ё || Ах; 





ЧЕ. 
Ед (Е) +Т, = = Х2бо | +Н, 60$ [22] Дх; (3) 


аЕ. 
Е, д (Е) +Т, и: = Х,с, 1+, со [©# || Ах, 





Для базовых динамических моделей, учитывающих деформационные смещения подсистем ин- 
струмента в пространстве и заготовки в плоскости, нормальной к оси её вращения, а также дина- 
мическую связь, объединяющую взаимодействующие подсистемы через процесс резания (2), (3), 
справедлива следующая система [3]: 








2х0 Хх 
т, —. +В, Г +сХ® =Е® м ); 
(4) 
2х () 
т и +В 7 + сх ® =Е® та х® }, 


где х® (Е) ={х, (Хх, (#),х, (у — деформационные смещения вершины инструмента в по- 
движной системе координат, привязанной к траектории движения суппорта; 
хе) (Е) =. (Е),х, (Е) — деформационные смещения заготовки в точке контакта с ней 


вершины инструмента В ПЛОСКОСТИ, нормальной к оси вращения шпинделя, 


т 0 0 п, № В, Ся Са Са 
т = 0 т 0, АВ=А, В, В, <= С, С,, 6;›| — параметры динамической 
0 о м, В, В, В, С С23 Сзз 


подсистемы инструмента (матрицы инерционных, диссипативных и упругих коэффициентов); 


т 0 
с — матрица упругости подсистемы инструмента, представленная в виде (1); т -| | 


во 
р -| бы матрицы инерционных и скоростных коэффициентов подсистемы заготовки; 
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И |^5, [6 + ДЕ, (Е -т}} ихи Хх Е (Х,) „Ед (Х, „Руд (х,} — вектор — функции динами- 
ческой связи, формируемой процессом резания. Очевидно, Е“) ={65, [6 +Л6, (#-т)] {их + 


ы 
+ д (Хз) „Е, д (№, . Здесь $, = соп$Ё — постоянное значение величины подачи суппорта 
в продольном направлении; &, = соп$Ё — постоянное значение глубины резания в предполо- 
жении, что вариации припуска в пределах оборота отсутствуют; ЛЁ, ({-Т) — периодическая 


по периоду вращения заготовки функция изменения припуска и (или) радиальных биений 
шпинделя; А, „ (Х,) — функция приращения сил, представленная в виде (2). 


Прежде всего дадим качественную характеристику системе (4). Ранее показано [2], что 
матрицы т,, В, с,, т, В и с являются симметричными и положительно определёнными. Бо- 
лее того, ориентация эллипсоидов жёсткости и диссипации в этих матрицах совпадают. Система 
(4) имеет стационарное решение, определяемое из условия Х; (Ё) = Х. (#-Т). Это установившие- 
ся вынужденные колебания относительно некоторой точки Х“®. Координата этой точки является 


решением системы (4) при правой части, равной Е® = р5„Ё, {Х/Х», Хз} = сОПЕ и 


Е) =р$>Ёь ен = соп$Ё . Это точка равновесия системы при отсутствующих вариациях при- 


пуска. Рассматриваемая система является линейной, для неё справедлив принцип суперпозиции. 
Поэтому в окрестности стационарной траектории формируются периодические движения 


Хх (Е =ХО(Е+Т), определяемые внешними заданными периодическими силами 


Е = 5, ЛЕ, (Е-Т ХХХ и Е =р5> ЛЕ, (#-Т){Х,Х»!. Обычно при определении влияния 


вариаций припуска на траекторию движений инструмента относительно заготовки ограничивают- 
ся рассмотрением этой составляющей. Однако, за счёт связи (2), входящей в правую часть систе- 
мы (4), стационарная траектория может потерять устойчивость. Для этого необходимо рассмот- 


реть уравнение в вариациях относительно траектории Х (9 (Е) =Хх@+х® (Е) . 


После замены переменных Х (ЕЁ) = хе (Е) х (ЕЁ) и для малых отклонений от стационарной 


траектории с учётом запаздывающих аргументов получаем уравнение в вариациях относительно 
стационарной траектории 





ах ах 
т; +И, -[+С,Х =0, (5) 
[611 [614 
где т,,№,,с, — соответственно суммарные матрицы инерционных, скоростных и упругих коэф- 


фициентов динамической системы резания, рассматриваемой в вариациях относительно ста- 
ционарной траектории, обусловленной вариациями припуска и изменениями жёсткости под- 
системы инструмента. 


,7 


Е 

и 
ооо 
аа 2:5 
20300 
озооэо 
зоооо 
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п № Х, Т,р5> п, 








[=М 
В, › - Х.Т.р5ь №» ПВ 5СоГ› | +.Н; с0$ | [(2/- 1) 0 ] > 
1=1 
1=М 
В, = Въ -ХЬТ.р5ь №3 №3 -Х. СОТ, | +У.н, 60$] [21-10 > 
1=1 
-Х, 7.05, 0 У | + у ы, с0$| (2/ - ее > 
1=1 
1=М 
—ХТ.р5» 0 -Х›соТ, | + Ур, с0$ [(2/ - 0: > 
№ [= 
> —Х, 7.05, 0 | 
<> -Х›Т.р5ь 0 
<> -Х-7.25$> 0 
<> В-Х,Т.р5 - СТ (1+ 60530), 0 
> —Х.Т.р5 й-СТ Е +с0$ [39 + 3" 
Г 1=М —_ | 
С: +Х.Р5р С: С;1+Х,60 | + У`Н, 0$ [(2/ -1) 9 > 
1=1 г. 
1=М д 
Хр бо | +» Н, 60$ [(2/ -1) 9 > 
1=1 = 
1=М у 
= отр. бо | +У`Н, 60$ [(2/ -1) 9 > 
1=1 = 
Х,05, 0 Х, С. | + у ы; С0$ [(2/ - а > 
1=1 
1=М 
Х.р5 0 У. 5. | +У.н, с | (2/ -99 > 
|. 1=1 
<> Х,25ь 0 й 
> Х.Р5» 0 
> Х:55$» 0 
<=> с(1+Нс0$30)+хХ,р5, 0 
> Х.Р5ь с 1 +Нс0$ [30 — 3" 


Если дополнительно учитывать вариации скорости подачи, например, обусловленные вза- 
имодействием суппорта с процессом резания, система (5) ещё усложняется. Пусть задана траек- 


тория скорости подачи в виде И, (Ё) =И„,+ДИ, (ЕЁ). Примем во внимание, что величина подачи 


=М 
ми , 0$ | (27-1) 9 ] > 


и" 
+ 














определяется интегральным п реобразован ием 


Е 


$ (Е) = [Ио +АМ, (2) 4 =5,+4$, (Е), 


ЕТ 
Е 


где 45, (#) = | ЛУ, (Е)4Ё — также Т_ — периодическая функция. 


ЕГ 
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Поэтому за счёт вариаций составляющей силы РЁ, имеет место изменение скорости пода- 


чи. Причём, эти изменения имеют также период Т’. Поэтому в системе (5) образуются дополни- 
тельные составляющие в матрицах скоростных коэффициентов и упругости, периодические с пе- 
риодом Т. Наконец, имеют место периодические составляющие вариаций скорости подачи, не 
связанные с периодическими изменениями величины припуска. Они обусловлены, во-первых, 
особенностями взаимодействия суппорта с направляющими станка через узел трения. Во-вторых, 
они определяются периодически изменяющимися погрешностями изготовления зацеплений в ме- 
ханической части привода, то есть всеми факторами, влияющими на точность станка. 

Важно подчеркнуть, что уравнение в вариациях относительно стационарной траекто- 

рии (5) имеет периодические коэффициенты, влияющие на её устойчивость, то есть они характе- 
ризуют не рассматриваемый ранее в динамике станков механизм потери устойчивости, связанный 
с параметрическим возбуждением. 
Анализ устойчивости системы резания с учётом формирования периодически изменя- 
ющихся коэффициентов. Если в системе (5) отсутствуют вариации припуска и жёсткость под- 
системы заготовки по её периоду вращения остаётся неизменной, то мы получаем хорошо иссле- 
дованный случай [3]. Потеря устойчивости в этом случае связана с формированием циркуляцион- 
ных сил и преобразованием симметричной составляющей матрицы скоростных коэффициентов из 
положительно определённой в отрицательно определённую. Здесь сосредоточим своё внимание 
на влиянии периодически изменяющихся коэффициентов на устойчивость системы. Для этого 
необходимо воспользоваться теорией Флоке, включающей определение матрицы монодромии и 
анализ мультипликаторов системы (5), преобразованной к форме Коши. 

Вначале положим, что периодические вариации припуска отсутствуют и жёсткость подси- 
стемы инструмента на порядок превышает жёсткость подсистемы заготовки. Тогда система (5) 
преобразуется в две скалярных подсистемы вида: 


2 
тбХ нь с (1+ ряпох = су; 
СЕ [919 (6) 
ау 
т =х, 
О”. 


где с, — жёсткость процесса резания; Г — постоянная времени стружкообразования. 

Уравнение (6) при Т =0 есть уравнение Матьё — Хилла. В отличие от этого уравнения 
система при и=0 может потерять устойчивость равновесия в зависимости от соотношения пара- 
метров Т и с» (рис. 2). 

Рассмотрим случай, когда при этих параметрах, но при и = 0, система является асимпто- 
тически устойчивой. Как и при исследовании уравнения Матьё — Хилла, введём в рассмотрение 
дополнительный безразмерный параметр п = /, , где в, = (с + С )/т — собственная частота 
консервативной системы при и=0. На рис. 3 приведён пример изменения областей устойчивости 
В ПЛОСКОСТИ (мп), позволяющий определить в зависимости от и критические диапазоны частот 


вращения шпинделя, при которых система теряет устойчивость за счёт параметрического воз- 
буждения. Области неустойчивости на приведённой иллюстрации отмечены точечными кривыми. 
Приведённые три диаграммы отличаются значениями постоянной времени стружкообразования 
Т , которые показаны на рис. 2 звёздами. Если и =0, то при всех рассматриваемых параметрах 


система является устойчивой. Однако по мере приближения параметров к границе области устой- 
чивости, области параметрического возбуждения системы расширяются и формируются новые 
лепестки неустойчивости в области низких частот вращения шпинделя. Таким образом, области 
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неустойчивости системы за счёт параметрического возбуждения зависят не только от частоты 
вращении шпинделя и модуляции жёсткости, определяемой коэффициентом п, но и от парамет- 


ров динамической характеристики процесса обработки. 
с, КГ/ММ 


1000 





| | | Область 
800 | | | неустойчивости 
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Рис. 2. Область устойчивости в плоскости параметров [Се | при н=0 


По данным, приведённым в работах [4, 5], частоты первых форм колебаний шпиндельных 
узлов токарных станков лежат в диапазоне 120—180 Гц. Например, при собственной частоте пер- 


вой формы колебаний шпинделя, равной @&, =120 Гц, на частоте вращения шпинделя 
2400 об/мин система потеряет устойчивость за счёт параметрического возбуждения по мере уве- 


личения 7, соответственно при ни =0,5, н=0,38 ин=0,24. 


п п 
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Рис. 3. Область параметрического возбуждения системы в 


плоскости параметров (м, П = @/&,) : а— с, =200 , кГ/мм, 








б } Т-0, < 6— с, =200, кГ/мм, Т =0,4.10*, ©; с— 





с, =200 , кГ/мм, Т =0,8.10“*,с 
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Если учитывать динамические свойства подсистемы инструмента и принимать во внима- 
ние только первую гармоническую составляющую вариаций припуска, то области неустойчивости 
расширяются ещё в большей степени (рис. 4). Здесь м,, и, — коэффициенты параметрической 
модуляции соответственно первой и третьей гармонических составляющих в соответствии с си- 
стемой (5). Матрицы коэффициентов подсистем инструмента и заготовки приведены в табл. 1 и 2 
соответственно. Параметры динамической характеристики процесса резания: с, =200, кГ/мм; 


р=300, кГ/мм? ; х= ле = (0,50, 0, 70, 0,51}. Частота первой формы колебаний, соответ- 


ствующая подсистеме заготовки, равна 120 с`". 


Й = 6/6, П= 6/6 
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Рис. 4. Области параметрического возбуждения системы при учёте деформационных смещений инструмента и 
вариаций припуска 








Таблица 1 
Исходные матрицы динамической системы резания инструмента 
т.103, кГс?/мм В, кГс/мм с ‚кГ/мм 
0,079 0 0 0,0159 0,0032 0,0021 1000 200 120 
0 0,079 0 0,0032 0,0225 0,011 200 2100 310 
0 0 0,079 0,0021 0,011 0,032 120 310 2500 

















За счёт взаимного влияния отдельных форм колебаний на различных частотах, во-первых, 
область параметрического возбуждения колебаний и потери устойчивости стационарной траекто- 
рии существенно расширяется, во-вторых, она смещается в сторону низких частот. 
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Таблица 2 
Исходные матрицы динамической системы резания заготовки 
т, кГс?/мм В, кГс/мм с, КГ/мм 
0,0414 0 0,0159 0,0032 400 0,0 
0 0,041 0,0032 0,0225 0,0 400 

















Заключение. Периодические изменения параметров, как динамической характеристики процес- 
са резания, за счёт вариации припуска или за счёт периодического изменения жёсткости подси- 
стемы обрабатываемой заготовки вызывает на определённых частотах вращении шпинделя поте- 
рю устойчивости процесса резания. Причём, области потери устойчивости изменяются при варьи- 
ровании параметров динамической связи, формируемой процессом резания. Этот механизм поте- 
ри устойчивости до настоящего времени в рассматриваемой предметной области не принимался 
во внимание. Кроме этого, необходимо отметить, что при изучении кинематических возмущений, 
действующих на систему резания, необходимо учитывать, что эти возмущения всегда вызывают 
изменение параметров системы. Поэтому изучение параметрических явлений позволяет объяс- 
нить ряд известных экспериментально фактов, например, увеличение вибраций инструмента от- 
носительно заготовки в вариациях относительно точки равновесия после возрастания частоты 
вращения шпинделя выше некоторого критического значения. Приведённые результаты имеют 
принципиальное значение особенно в тех случаях, когда рассматривается скоростная обработка 
резанием. Приведённые результаты также показывают, что уменьшение точности станка, связан- 
ное, например, с радиальными биениями шпинделя, расширяет область неустойчивого резания. В 
существующих представлениях о динамике процесса резания кинематические и другие погрешно- 
сти станка не влияют на устойчивость, и область устойчивости всегда расширяется по мере уве- 
личения скорости резания, что противоречит известным экспериментальным результатам. 
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